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We report a study of D*j(2317) + and D s j(2460) + meson production in B decays. We observe the 
decays B + -> D^j )+ D M0 and B° -» D^* )+ D ( * )_ with the subsequent decays D* J (2317) + -> D^tt , 
D sJ (2460) + -» £>+7, and L> sJ (2460) + -> D* + tt . Based on a data sample of 122.1 million BB pairs 
collected with the BABAR detector at the PEP-II B factory, we obtain branching fractions for these 
modes, including the previously unseen decays B — > D[*J + D* . In addition, we perform an angular 
analysis of D 3 j(2460) + — > Dfy decays to test the different U s j(2460) + spin hypotheses. 

PACS numbers: 13.25.Hw, 13.25.Ft, 14.40.Lb 



The unexpected observation of a narrow Dfir reso- 
nance with a mass of 2317MeV/c 2 was recently reported 
by the BABAR collaboration [l( and confirmed by the 
CLEO experiment 0. CLEO observed a second D* + tt° 
resonance with a mass close to 2460 MeV/c 2 Q , previ- 
ously suggested and later confirmed yj by BABAR. 
The Belle collaboration confirmed both resonances and 
found two additional decay modes for the higher-mass 
state -D+7 and Dfir + n~. These resonances are usu- 
ally interpreted as P-wave cs quark states I S El El El, 
although other interpretations [ljj, [ill [lj, [l3|, [lj] cannot 
be ruled out, and will be referred to in the following as 
2^(2317)+ and D s , 7 (2460)+ mesons. 

The new states were first observed in 
lisions. Their observation in exclusive B 



— » cc col- 

(*)+ 



decays allows additional properties of the D s j states 
to be studied: the D s ,/(2460) + -> L>+ 7 helicity angle 
distribution in B decays can be used to obtain informa- 
tion on the Z? s ,/(2460) + spin J [l^, and the measurement 
of the different branching fractions can help clarify the 
nature of these states. 

)(*)+ 

}(*)+£)( 



In this Letter we consider the D^j + production modes 

with the sub- 



B+ -> D { S *J + D^° and B° -> D sJ 
sequent decays L»* J (2317)+ -> L>+tt°, D sJ (2460) h 
£>*+7r°, and D sJ (2460)+ -> D+-/. Our intention is to ob- 
serve previously unseen decay chains, measure branching 
fractions for all channels, and determine the D s j(2460) + 
spin by means of an angular analysis. Charge-conjugate 
reactions are assumed throughout this paper. 

The measurements reported here use 113 Ad of data, 
corresponding to (122.1 ± 1.3) x 10 6 BB pairs, collected 
at the T(4S) resonance with the BABAR detector [l6j at 
the PEP-II asymmetric-energy B factory. 

We reconstruct D and Df mesons in the modes D° — > 



and 



K + ir , K + tt 7T°, K + ir ir + ir : D — > K + ir ir 
D+ -> 07T+ (<f> -> K+K-), K*_^K+ (K*° -> K-7T+). 
The reconstructed mass of the D and candidates is 
required to be within 2.5a (3a for K + n~ , K + 7r~Tr~ , and 
</>7r + ) of the nominal _D masses, where the D mass res- 
olution a, found in the data, is close to 12MeV/c 2 for 
D — > K + ir~ir° decays and varies from 5.3 to 6.3MeV/c 2 
for the other decay modes. 

The D* candidates are reconstructed in the decay 
modes D*+ -> D q tt+, D*° -► D°ir°, £>° 7 , and £>*+ -> 
£>+7. The mass difference Am between the D* and 
Z? candidates is required to be within 2MeV/c 2 of its 



nominal value Ha for D*+ -> D tt+ and D*° -> L>°7r° 
(10MeV/c 2 for D*° -> L>° 7 and D*+ -> L>+ 7 ), corre- 
sponding to about 5(JAm for £>*+ and 2dA m for Z3*° and 
D*+. 

(*)-|- (*) 

The selected pairs of D s and D candidates are 
combined with a photon or a ir° to form B candi- 
dates. The photon energy is required to be greater 
than 100 MeV. The neutral pions are built from pairs 
of photons with energies above 30 MeV and an invariant 
mass between 115 and 150 MeV/c 2 . A mass-constrained 
kinematic fit is applied to all the intermediate particles. 
In order to suppress combinatorial background, we re- 

quire the D invariant mass to be greater than 

2.3GeV/c 2 and 2.4GeV/c 2 for D and D* final states, re- 
spectively. Events compatible with being two-body de- 

cays B —y Dg D are rejected. 

We define a B signal region in terms of the beam energy 
substituted mass, toes = \/ s /4 — Pg 2 , and the difference 
between the reconstructed energy of the B candidate and 
the beam energy, AE = Eg — y/s/2, where s/s is the total 
energy in the T(AS) center-of-mass frame and Eg (p* B ) is 
the energy (momentum) of the B candidate in the same 
frame. We require 5.272 < toes < 5.288 GeV/c 2 and 
\AE\ < 32(40) MeV for tt° (7) final states. The width 
of the signal box is approximately ±3<r in toes and ±2tr 
in AE. We also define in the Ds + 7r°/7 mass spectra a 
signal region \m(D s + n° /^y) — m(D < s *J + )\ < 2.5er and a 
sideband region from 4ct to 12ct away from the nominal 
value, with m(D { *] + ) = 2.317 GeV/c 2 (2.460 GeV/c 2 ) for 
D+ tt° (D*+ 7r°, D+ 7). The resolution a = 8MeV/c 2 
(12 MeV/ c 2 ) for tt° (7) final states is obtained from sim- 
ulated signal events. 

The AE, toes, an d D^ + 7r° or Df~/ mass spectra of 

the selected events are shown in Fig. ^ for each of the 

(*)+ 

three D s j final states after combining the charged and 
neutral B — ► D^*J + D^ modes and summing over all the 

and D s *^ + decays. Data points in each plot show the 
distribution of one variable in the signal regions of the 
other two. We also show (cross-hatched histograms) the 
AE and toes spectra of events in the D[*J + sidebands. 
Only one B signal candidate per event, based on the 

smallest |A_E|, is entered in the Di*^ + ir and -D+7 mass 

(*)+ 

spectra and kept for further analysis. The D\j yields, 
masses, and resolutions, obtained from fitting a Gaussian 
signal function and an exponential background to these 
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FIG. 1: AE (left), m E s (center), and m(Z^* )+ ) (right) spectra for the B -> D ( *J + D w candidates: (a) D*j(2317)+ -> D+vr , 
(b) -Dsj (2460) + — > Ds + 7T°, and (c) D s j (2460) + — > Df'y. Data points in each plot show the distribution of one variable in the 
signal regions of the other two, as defined in the text. For AE and rriBS, the cross-hatched histograms are from events from 
the m(D[*J + ) sideband regions defined in the text. For the m(D[*J + ) plots, only one B signal candidate per event has been 
selected. Curves correspond to fit results. 



TABLE I: Event yields, reconstructed D *j masses and res- 
olutions in each final state for B — > D^*J + decays. 



Decay mode 


Yield 




Cm 






[MeV/c 2 ] 


[MeV/c 2 ] 


D* sJ (2317) + D M [DfTv ] 


88 ±17 


2317.2 ± 1.3 


5.9 ± 1.4 


D aJ (2460)+ £> W [D* s + ty°] 


112 ± 14 


2458.9 ± 1.5 


10.8 ±1.3 


D s j{2A&Q) + D { * ] [D s + 7 ] 


139 ± 17 2461.1 ±1.6 


12.1 ± 1.6 



spectra, are given in Tabled The measured resolutions 

are compatible with expectations from the simulation, 

(*)+ 

assuming zero intrinsic width for D s j . We have also 
confirmed that the yields obtained from fits to the m-ES 
and AE spectra are in good agreement with the yields 
fitted from the D s j mass spectra. 

The branching fraction measurement is based on the 
individual Dfn a , D* + tt°, and -D^7 mass spectra for each 

of the twelve D s D 71- /7 final states. As shown in 
Fig. |2 signals for B -> D [ *] + D^ are observed in all 
channels. The results of likelihood fits to these distribu- 
tions, using a Gaussian signal and an exponential back- 
ground function, are overlaid. In these fits, the Gaussian 
mean value is fixed to 2317 (2460) MeV/c 2 for L>* 7 (2317)+ 
(L> sJ (2460)+). Its width is fixed to 8 (12) MeV/c 2 for tt 
(7) final states, as obtained from the simulation and con- 
firmed with the data (Table The D[*J + event yields 



and the statistical significances are listed in TablelTTl The 
significance is defined as y— 2 ln(Co / C rnax ) , where C max 
and Co are the likelihood values with the nominal and 
with zero signal yield, respectively. A significance larger 
than 4 is observed for 10 of the 12 modes. 
(*)+ 

From the D 8 j event yields in the data, we compute 
cross-feed-corrected branching fractions, using the signal 
efficiency and the relative contributions from cross-feed 
between the different D S J decay modes as obtained 
from simulated signal events. The resulting branching- 
fractions are given in Table ILT1 together with the efficien- 
cies, including the intermediate branching fractions, and 
the internal cross-feed contributions. 

The dominant systematic errors come from the track- 
ing efficiency (1.3% per track), 7 and tt° efficiencies (2.5% 
per 7), the Am requirement on the D* + and D*° selec- 
tions (« 5% per D*), efficiency of the AE requirement 
(w 6%), Dl *J + mass resolutions assumed in the fits (5 to 
10%), and background fitting model (5%). We assume 
equal production rates for B + B~ and B°B° pairs and do 
not include a systematic error related to this assumption. 
The errors from the individual Z?W and Dl*^ branching 
fractions, as taken from |l7j . are given separately (Ta- 
ble[nj. They are dominated by the 25% relative error on 

B{D+ -» 07T+). 

From the measured branching fractions for B — > 
J D+ J (2460)!D ( * ) in the D*+tt° and in the L>+ 7 final states, 
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FIG. 2: £>+tt° (left), D s * + tt (center), and £+7 (right) mass spectra of the selected B signal candidates for the 12 D *j 
final states. Curves are the results of the fits described in the text. 



TABLE II: Event yields (including internal cross-feed 
ciencies (including intermediate branching fractions), 
error on B is statistical, the second is systematic, and 



contributions), number of events attributed to internal cross-feed, effi- 
and final branching fractions, B, for B — > D (*J + D (t) decays. The first 
the third is from the D and Df branching fractions. 





B mode 




Yield 


Cross-feed Efficiency (10 4 ] 


S(10" a ) Sif 


mificance 


B° 


^D* sJ (2317) + D- 


[D+A 


34.7 ± 


8.0 


0.3 


1.6 


1.8±0.4±0.3±o;4 


5.5 


B° 


^D* sJ (23l7)+D" 


[DtA 


23.5 ± 


6.1 


0.0 


1.3 


1.5 =b 0.4 ± 0.2±g;| 


5.2 


B+ 


-^D*j(23l7) + D° 


[DtA 


32.7 ±10.8 


0.3 


2.6 


1.0 ± 0.3 ± 0.1±g;a 


3.1 


B+ 


-> D* sJ (2317) 4 D*° 


[DtA 


17.6 ± 


6.8 


7.2 


1.0 


0.9 ± 0.6 ± 0.2±g;| 


2.5 


B° 


^D sJ (2460)+D- 


[d; + A 


17.4 ± 


5.1 


0.1 


0.5 


2.8 ± 0.8 ± 0.5lo;e 


4.2 


B° 


^D s j(2460)+D*- 


[d* s + A 


26.5 ± 


5.7 


0.0 


0.4 


5.5 ± 1.2 ± 1.0±i | 


7.4 


B+ 


^D sJ (2460)+ D° 


[D* s + A 


29.0 ± 


6.8 


2.2 


0.8 


2.7±0.7±0.5 + Jo;6 


5.1 


B+ 


^D s j(2460)+D*° 


[D* s + A 


30.5 ± 


6.4 


2.5 


0.3 


7.6 ± 1.7 ± l.StH 


7.7 


B l > 


~^D s j(2460)+D- 


[Dt-y] 


24.8 ± 


6.5 


0.5 


2.6 


0.8 ±0.2 ± 0.1±g;| 


5.0 


B° 


^D sJ (2460)+ J D*- 


[Dtj] 


53.0 ± 


7.8 


0.1 


1.9 


2.3 ± 0.3 ± 0.3±g;| 


11.7 


B+ 


^D sJ (2460) + D° 


[Dtj] 


31.9 ± 


9.0 


1.4 


4.1 


0.6±0.2±0.llo.i 


4.3 


B+ 


^D aJ (2460) 4 D*° 


[Dtj] 


34.6 ± 


7.6 


6.5 


1.7 


1.4 ± 0.4 ± O.3I0.3 


6.0 



we compute the ratio 

g(P sJ (2460)+ -> D+j) 
B(£> sJ (2460)+ -> D*+n°) 



0.274 ± 0.045 ± 0.020, 



where the first and second uncertainties are statistical 
and systematic respectively. This is compatible with the 
prediction from |fj( . 

We perform a hclicity analysis of the £> s j(2460) + state, 
using the decays B+ -> D aJ (2460) + lf and B° -> 
D sJ (2460)+D- , with D sJ (2460) + -► D+j. The helicity 
anglc 6h is defined as the angle between the D S J mo- 
mentum in the _B-meson rest frame and the D s momen- 
turn in the D s j rest frame. Since the Dj mass is cor- 



related with the helicity angle, the selection requirement 
m(D'y) > 2.3GeV/c 2 is omitted for the angular analysis. 
We perform m{D s n f) fits for five different cos(6h) regions, 
using the same fit functions and parameter values as in 
the corresponding branching fraction measurements. 

The resulting angular distribution, after applying cor- 
rections for detector acceptance and selection efficiency, 
is shown in Fig. [3] The predicted spectra for two dif- 
ferent assumptions for the D s j(2460) + spin, which have 
been normalized to the data, are overlaid. We exclude 
the J — 2 hypothesis (x 2 /n.d.f.=36.4/4) and find good 
agreement with J = 1 (x 2 / n - di -= 4 -0/4). A D sJ (2460) + 
spin J — is ruled out by parity and angular momentum 
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FIG. 3: Helicity distribution obtained from m{D s y) fits in 
bins of cos(#h) for data (points) in comparison with the ex- 
pectations for a D s j(2460) + spin J — 1 (solid line) and J = 2 
(dashed line), respectively, after normalizing the predicted 
spectra to the data. 

conservation in the decay Z? s j(2460) + — > Df-y. 

In summary, we have observed and measured the 
branching fractions for the decays B — ► £>*j(2317) + D^*^ 
(L>;, 7 (2317)+ -» L>+7T°) and B -» L> s j(2460) + _DM 
(UU/(2460)+ -» £>* + 7T°, D+7). The modes involv- 
ing a D* have been seen for the first time. The an- 
gular analysis of the decay B —> D s j (2460) + D with 
D Si7 (2460)+ -> L>+7 excludes J p = 2+ and supports 
the hypothesis that the D s ,/(2460) + is a J p = 1+ state. 
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